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In this paper, impedance characteristics are determined for discrete 
movement. It explores the changes in impedance characteristics of human 
arm during a complete task. This study considered 3D spatial movement for 
horizontal adduction and abduction. The human arm is considered as mass- 
spring-damper system and the modelling is done accordingly. The model is 


solved for 3 degree of freedom (3DoF) spatial movement usually used for 
daily work. Inertia, stiffness and damping factor are the impedance 
Keywords: characteristics considered in the model. Using position measuring device, 
the position data of elbow and wrist were obtained with respect to pre- 
defined references. The data were used to calculate velocity, acceleration 
: and force. Then, the impedance characteristics were determined by solving 
Mass-spring-damper the equation of motion of the mass-spring-damper system corresponding to 
Spatial movement different position of the wrist on the trajectory. These impedance factors 
Trajectory were then plotted to map the characteristics. The mappings were done for 
both total durations to complete task and for segments of the task. The 
impedance characteristics were observed to vary according to the movement 
of the arm. Significant changes of damping factor and stiffness were found 
to occur at the beginning and end of the arm movement. 
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1. INTRODUCTION 

Human arm movement is an emerging field which seeks to study how human arm works. In such 
research, dynamic impedance characteristics are mapped against time to have a picture of the dynamics of 
arm movement in space. The impedance characteristics are then used for designing robotic arm exoskeleton 
for rehabilitation. However, only few works have shown the dynamic changes of impedance characteristics 
for spatial movements. Humans have excellent abilities to manipulate objects as various task dynamics are 
adapted by the central nervous system. For example, for the first time due to some friction, one may have 
difficulty in opening a door. However, after one or two trials, the appropriate force is applied, and without 
difficulty one will open the door even without thinking about it. With the help of tools, many tasks performed 
by humans are essentially unstable [1]. Consequently, there is the need for continuous acquisition of 
additional skills since in unstable tasks dynamics, different initial conditions, neuromotor noise [2], or any 
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type of external perturbation can lead to unsuccessful and inconsistencies in the performance. Understanding 
how humans interact with the environment can provide insight into the neural mechanisms of adaptation and 
motor learning with potential application to computer animation of human motion. This will aid the 
development of robots collaborating with human operators and neuromotor rehabilitations. Conversely, 
humanlike control strategies are better employed in the improvement of robotic control especially with 
respect to safety [3]. 

Significant advances in medical robotics have been made recently (Healthcare robots 2021), and in 
micro and nanotechnology and teleportation [4], all of which rely on haptic interfaces to allow the human 
operator manipulate objects not directly accessible to the hand. The safety and performance of such systems 
significantly depend on motion stability i.e., the coupled stability of the arm interacting with the environment 
[5]. To properly control these types of systems, it would be useful to know the mechanical impedance. 
Mechanical impedance is the resistance to infinitesimal perturbations applied at the hand. The mechanical 
impedance of the human arm can be calculated from the restoring force to slight perturbations imposed in 
static positions [6] or during movement [7]-[11]. For a given stable or unstable interaction, it is possible to 
measure the impedance in the movements [12]. Though, this requires many movements, and in every 
dynamic interaction, it would be more useful to have a compact model to describe the force and impedance. 
Impedance is shown to depend on position force [6], [13] and instability [12] but no comprehensive model 
has been proposed so far. In order to learn more about human motor adaptation, works have investigated the 
adaptation to stable [14]-[16] and unstable [17]-[19] interactions produced by a haptic interface. 

Measurement of human arm impedance characteristics is an emerging field. This is due to numerous 
issues including the difficulty in the measurement of human kinematic data, the high complexity in the 
understanding of the working principle of human muscle, many hidden factors responsible for smooth arm 
movement. For a healthy normal human being, the arm movement needs to be very smooth, precise, and 
definite. It is obvious that for arm movement, not only the arm muscle is responsible, rather they are 
governed by the central nervous system (CNS) taking signal from the eyes and other sensory organs. When 
an arm of a human body is performing a task like picking an object or reaching to a specific point, it uses the 
five basic senses of the body. It uses the vision mainly to obtain information about the location of the object, 
and how it can be reached. The vision helps to gather enough information to perform desired tasks. The 
obtained information then computed by the brain and the whole arm system is given signal accordingly to 
perform the task. According to this information, the eye sends signals through the CNS to the arm muscles to 
pick or reach object accordingly. Most research have tried to reveal how the CNS controls the arm impedance 
and the changes in the impedance during arm movement. 

Disabled humans are unable to move their arm according to the CNS instruction due to damaged 
muscle cells. Due to different reasons, human beings could become partially disabled, and this has informed 
the development of different types of rehabilitation. This field combines the knowledge of biomedical 
engineering and mechatronics. Rehabilitation robots rehabilitate and help or assist the disabled people so they 
could perform any task without the help of any person. These robots are designed, fabricated and carefully 
attached to human arm so they could work as one with the human body. These robots help disabled persons 
to perform various functions in daily life. This way, disabled people would not have to rely on others to 
perform vital tasks. 

This paper presents the impedance characteristics of human arm movements and to map the 
impedance characteristics of human arm in spatial movements. This paper also reveals how impedance 
characteristics of human arm changes during a complete cycle of specific task. In this research, impedance 
characteristics of human arm considering spatial movement for horizontal adduction and abduction have been 
determined experimentally. This work is limited to horizontal adduction and abduction of human arm motion. 
The motion is 3D and involves fast movement as required for daily works. Simulations were done with 
different combinations to have better understanding of human arm impedance characteristics. To understand 
human arm impedance characteristics, it is also necessary to know different fields such as, biomechanics, 
biomedical engineering, robotics, and human anatomy. Analysis has been done by using MATLAB. These 
types of robotic systems are utilized in the calculation of human arm travelling trajectory. 


2. RESEARCH METHOD 

In this research, the movements (single and multi-degree of freedom) of human arm in space for 
various tasks are considered. The upper limb consists of three main parts namely, upper arm, forearm and 
hand. It is composed of three chain mechanisms, the shoulder girdle, the elbow and the wrist, whose 
association allow a wide range of combined motion. If the shoulder girdle moves so is the elbow and the 
hand is moved together accordingly. Here, both active and passive movement is investigated. The active 
movement refers to the movement when a person moves their arm by themselves while in passive movement, 
a person cannot move their arm as per their will except with some help. Suitable manipulators (single and 
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multi-degree of freedom) are constructed to move the arm and to measure the trajectory of the arm 
movements. The movements are executed as a compromise of the simultaneous, smooth interpolation of joint 
angles, discomfort minimization, and straight hand paths. Individual muscle and interaction between them at 
the time of movement is considered. This is based on various minima models which are applied to the on 
human arm movements. For example, minimum jerk motion, minimum torque change model, and minimum 
energy cost model. 

The experimental data (position, velocity, acceleration, and force) is analysed by using different 
variants of Recursive least-square method for mapping the impedance parameter. Finally, experiments are 
conducted to verify the proposed model. The steps taken in this research to achieve the stated research 
objectives are summarized thus: The research begins with comprehensive literature review to discover and 
explore the wide range of biomechanics, mechatronics and human anatomy, and consequently, analyse the 
related problem in arm impedance characteristics in spatial movement. The data of the elbow and the wrist 
are taken with the help of sensors. They are then analysed with MATLAB. The data of the active and passive 
movement are taken mathematically and plotted into graph for better understanding. The obtained plot is 
used for comparing all the proposed theories. The data are analysed by using different variants of recursive 
least-squares methods for mapping the impedance parameter. 


2.1. Experimental setup 

The Cartesian positional data in 3D space is measured using an array of fiber optic curvature sensor 
known as a Shape Tape. This tape is generally available with different numbers of sensors of different 
sensing zone lengths and with different leads and tip lengths. The schematic of a typical shape tape is as 
shown in Figure 1. The tape used in this experiment is 1810 mm, leads length is 1080 mm, the length of the 
sensing zone is 720 mm, and the tip is 10 mm. Its sampling rate is 110 Hz. In Figure 2, the shape tape device 
is attached to a subject. Data connector is connected to the device by the COM port and to the PC by USB 
connection. The arm tape used here is Shape Tape plus. For this case, the arm’s homing position was the 
palms facing down. Before acquiring the data, the shape tape reading was verified by moving the device in 
known environment. Also, it was checked to determine how many data points can be obtained. To be more 
specific, the length of the sensor portion is 760 mm. It gives the positional data with respect to the settings. 


Sensor pair intervals: 


ear 
L— Sensorzed— 
Zone 
Tip 


Figure 1. Schematic diagram of shape tape 


Figure 2. Shape tape device 


2.2. Experimental paradigm 

Six subjects including four females and two males participated in the experiments as shown in 
Table 1. Participants are ranged between 25-35 years old. None of the subjects had any sensory, perceptual, 
or any motor disorders. Subjects are not familiar regarding the research topic and had no prior experience of 
the experiment. Therefore, at first, they are given some instructions and asked to do the trials. In the 
experiment, both upper and lower arms are considered. The subjects were asked to move their hands in a 
parallel horizon with the surface of the table for some distance from it. Thus, the subjects were not dragging 
the mass, rather, they were carrying the weight. For uniformity of the arm trajectory, some indications were 
plotted on the travelling plane and the subjects were asked to move their arm according to the path. The 
parameter of subjects is given in the Table 1. 
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Table 1. Description of the subjects 


Subject Age Sex Length of arm, mm Length of lower arm, mm 
Sub 1 30 Female 60 33 

Sub 2 32 Female 54 30 

Sub 3 27 Male 59 32 

Sub 4 30 Female 61 36 

Sub 5 31 Female 56 30 

Sub 6 35 Male 56 32 


2.3. Data collection 

In this experiment, the subjects were asked to sit on a fixed chair to avoid unwanted body 
movements as shown in Figure 3. The middle point of the arm was fixed on the table. According to the 
paradigm, the subject shown in Figure 4 started from point | to point 2, and again, returned back to point 1. 
This whole procedure is considered as one full cycle. This one cycle was repeated for three times, but not 
continuously. Due to muscle fatigue and stress, the arm was given some space to take rest. The whole 
procedure was done periodically for three days. The considered movement is an active daily routine task with 
some external force applied by handling an object of 1.5 kg weight. 


Figure 4. Experiment pathway 


2.4. Data validation 

For the validation of the measured data, at first a kinematic model based on the inverse kinematics 
for the arm movement on the horizontal surface is developed and compared with the measured data. 
In Figure 5, the black dotted areas are the coverage of the arm where the arm can swing. The dark line at the 
right bottom corner of Figure 5 is the measured Cartesian data of the wrist. The arms limited movement is 
defined by considering the human arm movement limitation. The coloured line indicates the actual positional 
change of the wrist position. The sample is indicated by the colour bar of the right side which changes colour 
from blue to red. From the figure, it is shown that the measured data exist in the area where the arm can 
move. Also, the angular velocity at the shoulder and elbow are also estimated and compared with the 
theoretical results obtained from inverse kinematics [20]. 
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Measured data were used to calculate the angular position of shoulder joint and elbow joint using 
inverse kinematics [21], [22]. Theoretical angular position of the shoulder and elbow were compared with the 
calculated data. Data with minimum error were used for further calculation. A sample comparison is shown 
in Figure 6. The angular change in the shoulder is shown in above and the lower one is for angular change in 
elbow. In Figure 6, the black line indicates the actual angular position change, the blue line is the angular 
positional change computed using inverse kinematics and the red line is the error calculated by least square 
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Figure 6. Comparison of angular changes between actual measured angle and the computed angles 


2.5. Mathematical model 
The applied mathematical model in this study is the same as the dynamic equation of motion for 


mass-spring-damper system. To acquire the kinematic data, Measured Inc., Shape Tape device was used. 
This device gives the positional value of the desired Cartesian coordinates against time. The developed arm 
model for the measurement of the dynamic impedance is as shown in Figure 7. The hand impedance model 
can be expressed for the end-point level when the human arm is under a stable posture by (1). 


MX + BX + KX = F (1) 
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where M, B, KE R°* represent the hand inertia, viscosity and stiffness matrices respectively. 


Upper arm 


Lower arm 


Figure 7. Arm model for dynamic impedance measurement 


Xe R? and Fe R° are the end point position and force vector, respectively. The movement 
considered in this study is fast movement where some of the terms can be compensated [23]. Here the forces 
are calculated from the Newton’s second law. The position of the hand as well as the end position is 
measured by Measured Shape Tape® device. Therefore, as shown in (1) is transformed to an identifiable 
form, with respect to the impedance matrices M, B and K as given in (2). 


YI =F (2) 


where, IE R” is the parameter vector to be identified, which is given by (3) while the impedance matrix 
functions are given in (4) to (6). 


1 = [M] BT KT M? BT KT MT BT KT] (3) 

M = [M] M] MT] (4) 
—[ĪpRTpTpT]T 

B = | B] By B7] (5) 

K = [KT KUKI] (6) 


Y in (2) is a (3x25) matrix defined by 


LT Oix9 Oixo 
Y = | 0x9 LT 01x9 
Oix9 01x9 LT 


where, L € R? is defined by 
L= [xe XT XTT 


Oiis a zero matrix of size (1x9) X, Ë vectors representing end-point velocity and acceleration. 
According to (2) is used for computing I by using least-square method and Y is calculated by arm 
end-point positional data of one complete cycle. 


3. RESULTS AND DISCUSSION 

The main intention of this paper is to investigate the impedance characteristics of the human arm for 
designing exoskeleton for mimicking the human arm movement. In this study, the movement considered is 
3D and the duration of the movement is about 2 s. As the movement is not discrete point-to-point movement, 
it is not convenient to represent the impedance characteristics by ellipses. That is why the impedance 
characteristics against the time (sample) are plotted to see the dynamic changes during the whole process. 
In the plot, the impedance changes during the complete three cycles are presented. Experiment was 
conducted using six subjects. Due to the different ergonomic property, they have completed the assigned task 
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in different time and covered different workspace. This can be visualized in Figure 8. The data obtained from 
the experimentations are plotted against both with respect to sample and normalized time to see how the 
impedances differ from subject to subject and how the impedances change with the arm position with respect 
to time. A sample of wrist positions of six subjects are plotted against time in Figure 8. From this figure, it is 
seen that different subjects perform the same task in different time duration. The wrist positions change is 
plotted in normalized time while it is observed that all the subjects moved their arms following similar trend. 
It is also seen that subject 6 moved his arm lesser than the others, where subject 2 shows maximum 
movement. 

Figure 9 shows the three cycles of position X for six subjects. Downward curve represent the wrist 
movement in the forward direction while during the upward curve represent the movement of the wrist in the 
reverse direction. For all the subjects, the trend of X positions is similar for both forward and reverse 
motions. It is also seen that the X position maintains a smooth curve all the time which implies there was no 
jerk during this cyclic motion. The reach of the X position was different for different subjects due to their 
arm lengths and other argonomics data. 
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Figure 8. Position change along X direction among the subjects 


Position, mm 


Figure 9. Position plot against normalized time 


3.1. Inertia componenets 

The calculated inertia componemts are shown in Figure 10. From the figure it is seen that inertia 
change occurs during starting position, turning position and ending position. Moreover, considerable amount 
of change is seen in the X and Y diagonal components of inertia. As the Z directional movement is very 
small, the inertia change for all the Z diagonal components are very small except at starting position. 
Generally, for a solid object the inertia remains constant during the motion, but in case of muscle with 
respect to muscle movement muscle shape will be changing. Due to this reason, it changes its center of mass 
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and hence in every arm moment their inertia changes. It is shown that the trend is zigzag, and it ranges from 
20 N.s?/m to -20 N.s?/m, though around normalized time of 0.9, it gets out of 40 N.s?/m and -50 N.s?/m. In 
Figure 10 the considerable amount of change in X-diagonal component of inertia is observed at the ending 
position i.e. starting from position 1 to position 2 and again tending to return to position 1. From Figures 10, 
the maximum inertia is observed by subject 4 and minimum is by subject 1. The magnitude of the inertia of 
X diagonal component is larger than the Y and Z diagonal component of inertia. 
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Figure 10. Inertia componemts of all subjects 


3.2. Stiffness componenets 

The stiffness components of all dimension are shown in Figure 11. The maximum stiffness change 
along X component is observed in position 4 when the subject changes its movement direction and found by 
comparing with the position graph in Figure 9. The maximum magnitude is 7 N/m, and the minimum 
is -7.5 N/m while the maximum deviation is observed by subject 4 (Top graph). The maximum stiffness 
change is observed by subject 3 while the minimum by subject 1. The maximum stiffness change along Y 
component is less than X component. The maximum magnitude is 6 N/m, and the minimum is -4.5 N/m 
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while the maximum deviation is observed by subject 6. The maximum stiffness change is observed by 
subject 6 while the minimum also by subject 1. The stiffness changes along Z component are less than X and 
Y components as less movements in Z direction. 
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3.3. Damping coefficient componenets 
The damping coefficient components of all dimension are shown in Figure 12. The magnitude of Y 
diagonal component of damping factor is much smaller than the X diagonal component of the damping 
factor. The maximum and minimum points are observed by subject 5. For the Z diagonal damping factor 
observed during the starting and end of the motion, the average values lie between 100 N.s/m to 350 N.s/m 
for all subjects except subject 3 with the highest value of 500 N.s/m. For the minimum, it is not the same 
average value for all subjects except subject 3 with around -70 N.s/m to -350 N.s/m and highest value at 
-270 N.s/m. The average and standard deviation is around 20 N.s/m to 40 N.s/m though with the peak for 
subject 4 and subject 3 at 67.70 N.s/m and 71.17 N.s/m, respectively. 
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In this chapter the variation of impedance parameters is shown for different subjects with respect to 
time. From all the subjects it is seen that damping factor and stiffness are the mostly effected parameter 
during 3D hand movement. In all the studies it was observed that the magnitude changes of damping factor 
and stiffness are always high as compared to inertia. The magnitudes obtained for the impedances have 
shown similar trends as seen in other studies related to this field [24], [25]. 


A 
a A Subl 
> J + > 
eR Sub2 
:2 x  Sub3 
5A Sub4 
= 
> s! x Subs 
8 =F + * i "ow r " m + Subó A 
7 A " Mooca = n r 4 
eo : < + +x — + =| 
be Ar t+ Rah * * + * 
oS) 
av, = + 
a. 
5 al 
z 
A i - 
9 0.1 0.2 0.3 04 . 0.5 — 06 0.7 0.8 09 1 
Normalized Time 
ee Subi 
ga Sub2 
Z2 Sub3 
Rott # Sub4 P y 
~~ ae BA Í 2 # + Subs 
we OF A A + / < Â = A 
Sle. = sus = 1 + Subé 
Zot ss z~“ + #\ 4} FAEN SAIAS 
te i + * = a = A 4 
— 7 d TE Aki a 
So A A = à a ++ + = 
SF A 
s i + “a 
AS 
Q 
9 01 02 03 04 05 06 07 os 09 | 
Normalized Time 
A Sub! 
a. ^- Sub2 
zS Sub3 
P S = Sub4 
So ` + Subs 
Sa, =. Sub6 
=o $ 
vor 
o am # y # = 
A = + + + 
be ols Fee Se H t EOE + > re st). A . 
oe F 3 a RR k 
Ss -i + M + 
Amh y + + 
a x fe 
22 bi 
AN 
9 01 02 03 04 05 06 07 os 09 1 


Normalized Time 


Figure 12. Damping factor of all subjects 


4. CONCLUSION 

The main purpose of this paper is to map the human arm impedance characteristics during daily 
active movements. In this research, 3D human arm movement was considered. Human arm was considered as 
mass-spring-damper system. The model for mass-spring-damper system was used to map impedance 
characteristics. The position, velocity, acceleration, and force were used to solve the model for impedance 
parameters e.g. inertia, damping factor, and stiffness. These impedance parameters were plotted against time. 
In this paper, positional data were measured by Measuring’s Shape Tape device. During the experimentation, 
six subjects of age of 25 to 35 participated. Among them, four were female and two were male. Having 
analysed the impedance characteristics mapping, a mass-spring-damper system is found as appropriate model 
for mapping the human arm impedance characteristics. Mapping of dynamic impedance characteristics are 
done by plotting these characteristics against time. It is observed that the amounts of changes in damping 
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factor and stiffness are considerable whereas very little change is observed in inertia. Mainly, these changes 
occur due to the arm movement. 
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